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Techniques for studying the performance of Naflex pressure-assisted seals in the Ares-I 
Upper Stage liquid hydrogen tank manhole cover seal joint are explored. To assess the 
feasibility of using the identical seal design for the Upper Stage as was used for the Space 
Shuttle External Tank manhole covers, a preliminary seal deflection analysis using the 
ABAQUS commercial finite element software is employed. The ABAQUS analyses are 
performed using three-dimensional symmetric wedge finite element models. This analysis 
technique is validated by first modeling a heritage External Tank liquid hydrogen tank 
manhole cover joint and correlating the results to heritage test data. Once the technique is 
validated, the Upper Stage configuration is modeled. The Upper Stage analyses are 
performed at 1.4 times the expected pressure to comply with the Constellation Program 
factor of safety requirement on joint separation. Results from the analyses performed with 
the External Tank and Upper Stage models demonstrate the effects of several modeling 
assumptions on the seal deflection. The analyses for Upper Stage show that the integrity of 
the seal is successfully maintained. 


I. Introduction 

T HE objective of this paper is to perform analyses of the Ares-I US LH2 tank manhole cover seal joint to predict 
relative deflections between the cover and the LH2 tank dome-cover flange to assess the feasibility that the 
Naflex seal used with the ET manhole covers can successfully be used for the Ares-I design. Ares-I is the two-stage 
crew launch vehicle for NASA’s Constellation Program. Ares-I consists of a solid-fueled first stage, and a liquid- 
fueled Upper Stage (US) that usesliquid hydrogen (LH2) and liquid oxygen (LOX) to power asingleJ-2X engine. 

Theinsulated US structure includes two tanks to hold the LH2 and LOX at cryogenic temperatures. Thedesign 
of the domes includes “manholes” to provide access to the inside of the tanka During final assembly, covers will be 
placed over these manholes To prevent leaking of the liquid propellant or ullage gas when the tanks are filled, seals 
will be placed between the manholes and the covers 

In the Saturn and Space Shuttle Programs, Naflex seals were used to seal the joints between the manholes and 
the covers in the propel I ant tanka In this paper, a preliminary assessment is presented of the feasibility of using the 
Space Shuttle External Tank (ET) manhole cover Naflex seal design in the Ares-I US manhole cover joints As part 
of the assessment, a seal analysis technique different from that used by the ET Project is explored. 

Naflex seals are spacer-type, deflection-activated, pressure-assisted, static seals The seal configuration uses a 
cantilevered, deflection-loaded primary seal and a simple gasket-type secondary seal. A cross-sectional view of the 
seal is shown in Figure 1(b). The cantilevered beams or legs form a fork on the pressure side of the seal. In the 
uninstalled configuration, the tips of the cantilevered legs protrude beyond the rest of the seal in theX direction (see 
Figure 1(b)). The seal is installed in the joint between two flat flanges, and during joint assembly, when the bolts 
are preloaded, deflection of the legs provides the initial contact load to accomplish sealing at the primary seal-flange 
interface, creating the primary sealing surface. Pressurization of the joint forces the legs even tighter against the 
sealing surface(Anon, undated; Anon, ca 2000; Anon., 2001; and Robbins and Ludtke, 1964). In addition, the area 
around the bolt hole provides a secondary sealing surface as it is also placed in compression due to the bolt preload. 
Naflex seals are commonly fabricated of Inconel 718 with a Teflon enamel coating. The manhole cover seal design 
used by the ET Project has a heritage traceable to the Saturn S-l I stage. The Naflex seal under consideration in this 
paper i s shown i n Fi gure 1 . The seal has a 92-hole bolt ci rcle around the ci rcumf erence. 
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Figure 1. Naflex seal. 


The paper i s organi zed as f ol I ows. Fi rst, anal yti cal background of the axi symmetri c anal ysi s techni que and test- 
verification program employed by the ET Project is provided. Next, the analyses performed for this paper are 
described. Three-dimensional finite element modelsfor the ET LH2 tank manhole cover seal are correlated with the 
heritage ET seal test results Then, the modeling practices established by the ET analyses are used to construct 
modelsfor the Ares-I USLH2tank manhole cover, and the seal deflection analyses are presented. 


II. Analytical Background 

The sealing performance of theET Naflex seals was originally qualified by similarity (Anon., 1976), presumably 
to heritage Saturn-era seals Subsequent computer simulations of the manhole cover joints for both the LH2 and 
LOX tanks of the ET were performed during the 1980s (most recently, Pilet and Geiman, 1997) using BOSOR5 
(Bushnell, 1974 and 1976). BOSOR5 is afinite difference energy minimization code that performs large-deflection 
axi symmetric stress analysis, small -deflection nonsymmetric stress analysis, and buckling analysis For an applied 
pressure loading, relative deflections between specified nodes that represented the seal, the tank dome (the 
manhole), and the manhole cover were computed. These deflections were compared to the seal deflection limits 
levied in the Saturn-era seal design specification. This specification on seal deflection requires that the primary seal 
not leak unacceptably at established temperature and pressure conditions and with atotal joint axial deflection (joint 
opening) of Furthermore, this specification is usually interpreted as a one-sided seal deflection such that each 
flange of the joint opens no more than AspJ2. The worst-case deflection predicted by the ET analyses was at the ET 
LH2 aft manhole cover joint and exceeded Ajpec. A test program was then conducted to validate the BOSOR5 
modelsand also to determine the maximum seal deflection that can be present before sealing capability is lost. 

The test program verified the sealing capability of the seals for the ET tank configurations, temperatures, and 
pressures, and the BOSOR5 analysis models were correlated to the test results (Gillespie, 1988). The design 
specification for seal deflection was determined to be overly conservative, and a new seal deflection allowable for 
the ET joint configuration A CT was derived. 

All subsequent analyses of seal deflection under ET flight loads conducted by the ET Project continued to use 
the BOSOR5 modeling technique with the derived seal deflection allowable A ET as the success criterion. 
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III. Analysis 

Seal deflection analysis of the US LH2 tank manhole cover seal joint initially attempted to use the same BOSOR5 
axisymmetric modeling technique employed by the ET Project because the US LH2 tank manhole cover design is 
si mi lar to that of the ET, and the ET Naflex seal models had been anchored to test data 

However, after attempting this analysis, it became apparent that a different modeling technique was necessary 
because 

• The ET manholecover BOSOR5 models were created and run in an older, customized version of BOSOR5 
that is no longer supported nor compatible with newer versions. 

• Newer versions of BOSOR5 may not provide the most efficient method for this type of analysis due to the 
axi sy mmetri c si mpl ificati ons of the bolted j oi nt. 

Therefore, three-dimensional (3D), symmetric wedge finite element (FE) analyses are performed using the 
ABAQUS 3 commercial FE software. All of the analyses presented in this paper are nonlinear analyses. The bolt 
preload is applied during the first step of the analysis, and then the pressure load is applied during a second, 
subsequent step. 

The 3D FE analysis technique has several perceived benefits. 

• The LH2 tank dome and manhole cover designs of the ET and US are si mi lar but not exactly the same; a 
new FE model built uniquely for US allows the behavior of theLH2 tank manholecover joint to be better 
understood. 

• The design of the LOX tank cover on US is significantly different from ET, requiring new modelsto 
eventually be created. 

• A three-di mensi onal anal ysi s methodol ogy al I ows some of the assumpti ons associ ated wi th the 

axi sy mmetri c model i ng to be rel axed and the features of the seal j oi nt to be model ed expl i ci tly . Because 
more features of the seal joint can be modeled explicitly, more options are available for interrogating the 
anal ysi s resul ts and assessi ng the performance of the seal . 

A. 3D Finite Element Seal Joint Analyses 

As a first step in a methodical approach to the analysis, symmetric wedge models for the ET LH2 tank manhole 
cover seal joints are correlated with the ET seal test results (see Gillespie, 1988). The 3D FE modeling practices 
establ i shed wi th the ET model s are then used to construct an A res- 1 U S L H 2 tank model to assess the manhol e cover 
(MHC) seal performance under flight envi ronments. 

1. ET Seal Test-Analysis Correlation 

a. ET Seal Test 

The ET seal test was designed partly to demonstrate the ability of the BOSOR5 modelsto predict measured seal 
deflections Because the worst-case seal deflection predicted by the BOSOR5 analyses was at the LH2 tank aft 
manholecover joint, the ET seal test utilized a flight-ready LH2 tank aft MHC. A simplified view of the ET LH2 
tank aft MHC is shown in Figure 2. This MHC and a flight -quality Naflex seal were mounted to a rigid fixture tool 
as depi cted i n Fi gure 3. The test was conducted at room temperature. 

The test article was pressurized in 25-percent increments up to a specified target pressure. At each pressure 
increment, the vertical defection of the MHC was measured at five locations, D1, D2, D3, D4, and D5, as shown in 
Figure 4(a). Locations D1 , D3, and D5 are each in line with a bolt, while locations D2 and D4 are between bolts. 
The deflection was measured, as shown i n Figure 4(b), on the top surface of the M HC at the radi us correspondi ng to 
the seal primary sealing surface, flps& It was assumed that no relative motion occurred between the outer and inner 
surfaces of the MHC; the rigid fixture was assumed to not deform, and all of the vertical deflection was due to 
translation of the M HC. Therefore, the measured vertical deflection on the outside of the cover was assumed to be 
an accurate measure of the seal deflection. The finite element analyses performed for this paper validate this 
assumption, and hence, for conveniences presentation, the results will be presented as seal deflections. 

The seal deflections vs. pressure for the ET test are shown in Figure 5. In Figure 5, the normalized seal 
deflections at each of the locations D1-D5 and the average seal deflection for all of the locations are plotted. The 
seal def lections are normalized by the average seal deflection at the target pressure. 


3 ABAQUS is a registered trademark of Dassault Systemes. 
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(b) Cross-sectional view 


Figure 2. Simplified view of ET LH2 tank aft MHC. 
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Normalized Seal Deflection 




(b) Cross-sectional view: Radial location of 
deflection measurements 


(a) Top view: Gauge locations around theMHC 
(dra/vingfrom Gillespie, 1988) 

Figure 4. Locations of vertical deflection measurement during ET seal test. 



Figure 5. Seal deflections for the ET seal test. 
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b. 3D Finite Element Model 

A 3D solid FE model of the ET LH2 tank aft MHC iscreated by exploiting the near symmetry of the MHC. A 
wedge section of the M HC that encompasses a si ngle bolt, or 3.91 3°, as shown i n Figure 6(a), is used. A close-up of 
the symmetric wedge FE model is shown in Figure 6(b). The local region around the bolt hole is modeled with 3D 
solid elements, and the region away from the bolt hole is modeled with shell elements. The two regions are joined 
using a shell -to-sol id coupling constraint for smooth transition of loads across the shell/solid interface. This practice 
enables the use of shell theory, which is appropriate for the tank configurations, in conjunction with solid modeling, 
whi ch i s requi red to capture greater detai I at the seal j oi nt. 

Similarly, a 3.913° wedge FE model of the Naflex seal iscreated. Thismodel is shown in Figure 7. Theseal is 
modeled entirely with solid elements. The Teflon coating is not explicitly modeled; rather, the effects are modeled 
using a coefficient of friction in the seal contact definition, which isdiscussed later. 


I 



(a) Symmetric wedge of aft M HC 



(b) Close-up view of 3.913° wedge finite element model 
with local cylindrical coordinate system 

Figure 6. Symmetric wedge of ET LH2 tank aft MHC. 



Figure 7. Finite element model of 3.913° wedge of Naflex seal. 
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The FE model of the ET test set-up is shown in Figure 8. The model consists of three main parts: the MHC 
incorporating both shell and solid modeling, the seal, and a 3.91 3° wedge to represent the rigid fixture. Because the 
test fixture is assumed to be rigid, the bottom surface of the fixture model isfixed (i.e., u r = u„ = u z = 0). Symmetric 
boundary conditions (i.e., u e = 0) are applied to the 9 = 0 and 9 = 3.913° cut planes. The preload is applied to the 
bolt during the first analysis step. The pressure is applied along the inside surface of the MHC and around the 
pressure si de of the seal I egs, as shown i n Fi gure 8, duri ng the second anal ysi s step. Contact i s def i ned between the 
seal and the M H C and between the seal and the f i xture. M ore detai I s of the contact model i ng i n these two regi ons 
will be presented later. In addition, contact is defined between the lip of the MHC and the outer surface of the 
fixture (see Figure 8). 

The modeling of the bolted joint is achieved by the use of beam elements to represent the bolt and “spider” 
constraints to represent the interaction of the bolt with the connected parts An example of the bolted joints in the 
ET MHCs is shown in Figure 9(a). These joints consist of a bolt, a corresponding threaded insert, and a washer. 
The representation of the bolted joint in the FE models is shown in Figures 9(b-c). The bolt preload is applied using 
the ABAQUS bolt load capability: a pretension section is defined on the bolt, the preload is applied as a force to the 
pretension section, and then the bolt length is fixed. On the MHC side, the spider constraint couples the 
displacements of the chosen nodes on the outer surface, representing the size of the washer bearing surface (see 
Fi gure 9(b)) , to the moti on of a si ngl e node correspondi ng to the bolt head. On the fi xture si de, the spi der constrai nt 
couples the displacements of the chosen nodes on theinside surface of the threaded insert hole (see Figure 9(c)) toa 
single node corresponding to the bolt end, representing the bolt-engagement BC. 

The solid elements are the ABAQUS C3D8I elements, which are 8-node brick elements with additional internal 
freedoms introduced as incompatible modes for improved bending behavior 4 * * 7 . In general, in the MHC, there are 
eight brick elements for every half inch of thickness. The shell elements are the ABAQUS S4 elements, which are 
4-node quadrilateral elements, and the beam elements are the ABAQUS B31 elements, which are 2-node 
Timoshenko beam elements. Overall, the test set-up FE model contains 16208 elements: 8 beam elements, 488 shell 
elements, and 15712 brick elements. 

Linear, elastic material properties are assumed for all of the materials The materials used for the current 
analyses are provided in Table 1. To simulate the rigid test fixture, the modulus is arbitrarily set three orders of 
magnitude higher than the moduli for the other materials 



Pressu re 


Naflex 
Rigid Fixture 


Symmetry BC along cut planes: 
Solid Elements: u ( , = 0 
Shell Elements: 0 r = 0 Z = 0 


Fixed BC along bottom surface: 
u r =u e = u 2 = 0 


Bolt 


Contact between 
MHC and fixture 


Figure 8. Finite element model of test set-up. 


4 The C3D8I solid element has 8 nodes with three degrees of freedom per node plus 13 additional element variables 

associated with the incompatible deformation modes The estimated total number of unknowns for a given finite 

element mesh using C3D8I elements is roughly equal to three times the number of nodes plus 13 times the number 
of C3D8I solid elements 
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(a) Bolted joint in ET MHC 
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(b) Spider constraint on M HC side 


(c) Spider constraint 
on fixture side 


Figure 9. The bolted joint. 


Table 1 . Linear, elastic material propertiesfor the ET seal test set-up FE analyses. 


Part 

Material 

E, psi @ RT 

v @ RT 

MHC 

Aluminum 2219 

10.5x10 6 

0.33 

Seal 

Inconel 718 

30x1 0 B 

0.3 

Bolt 

A286 Stainless Steel 

30x10® 

0.29 

Fixture 

n/a 

10x10“ 

0 


The modeling of the contact between the seal and the M HC and between the seal and the fixture requires special 
handling. As discussed previously, the fixed boundary condition (which removes the rigid body modes from the 
problem) is applied to the bottom of the fixture, and the M HC is connected to the fixture via the bolt modeling and 
spider constraints However, the seal, which sits between the MHC and thefixture, isan unconnected region; i.e., it 
is not physically connected to either the MHC or the seal, and hence could cause numerical singularities. A cross- 
secti onal vi ew of the seal contact surfaces i s shown i n Fi gure 1 0. Inthefigure, initial contact at the pri mary seal i ng 
surface (PSS) and an initial gap at the secondary sealing surface (SSS) are shown. Note that the surface area of the 
SSS is much greater than the surface area of the PSS. When modeling contact in ABAQUS, any surfaces that are in 
initial contact with each other allow the corresponding nodal entries in the global stiffness matrix to be populated 
during the initial set-up (the zeroth solution increment), thus reducing the chance for numerical singularities to 
occur. I n order to take advantage of thi s, two separate contact i nteracti ons are def i ned for the PSS and the SSS. 

For both contact interactions, a finite-sliding, surface-to-surface formulation is used. The finite-sliding 
formulation utilizes the true representation (rather than a linear approximation) of the master surface for the slave 
surface i nteraction. The surface-to-surface formulation enforces contact i n an average sense over the slave surface, 
rather than enforcing contact at the individual nodes on the slave surface. Finite-sliding and surface-to-surface are 
the most general contact formulations available in ABAQUS, and they tend to yield better solutions for a wide range 
of problems than the other formulations (Anon, 2009). The constraint enforcement method in the direction normal 
to the contact surfaces is defined as “Hard” contact. In Hard contact, zero clearance between the contact surfaces is 
strictly enforced whenever contact pressure is detected. The constraint enforcement method in the direction 
tangential to the contact surfaces uses a penalty formulation. Because the seal is not physically connected to the 
other regions of the model, frictionless contact causes numerical difficulties, and a converged solution cannot be 
obtained. Therefore, a coefficient of friction is assigned for contact in the tangential direction. The use of a 
coefficient of friction, while not necessarily conservative, is not unreasonable in these problems because the Teflon 
coati ng on the seal provi des some resi stance to si i di ng. 

For the contact i nteracti ons between the seal and the M HC, the master surface i s def i ned as the bottom surface of 
the MHC, as highlighted by the green line in Figure 10. For the PSS, the slave surface is defined as the top surface 
of the seal ti p, as i ndi cated by the yel I ow I i ne i n Fi gure 1 0. To ensure that the two surfaces are recogni zed as bei ng 
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Figure 10. Seal contact surfaces. 


in initial contact, and hence populate the global stiffness matrix, the ABAQUS option to adjust the slave node initial 
posi ti ons i s used. For the SSS, the si ave surface i s def i ned as the top surface of the seal , as i ndi cated by the magenta 
line in Figure 10. To maintain the initial gap, the slave node initial positions are not adjusted. Identical contact 
i nteracti ons are def i ned between the seal and the f i xture. 

Initial analyses performed for this model did not converge. Because the surface-to-surface formulation enforces 
contact in an average sense over the slave surface, it generates unsymmetric terms in the global stiffness matrix. 
Therefore, the use of the unsymmetric matrix storage and solver was required to obtain converged solutions. 

c. Results 

The baseline FE model of the ET seal test set-up, hereafter referred to as the Baseline-ET model, is shown in 
Figure8 and is characterized by thefollowing: 

• The bolt preload i s the preload used i n the ET test. 

• The pressure I oad i s the target pressure from the ET test. 

• The ri ng of nodes for the spi der constrai nt on the M H C encompasses two ri ngs of el ements, as shown i n 
Fi gure 9(b) , thus model i ng a beari ng surface def i ned by the washer outer di ameter. 

• The nodes for the spi der constrai nt on the fi xture i ncl ude al I of the nodes on the bottom surface of the 
threaded insert hole, as shown in Figure 9(c). 

• The coeffi ci ent of f ri cti on for the contact between the seal and the M H C and between the seal and the 
fixture is (x = 0.1. 

• The contact between theMHC and the fixture (see Figure 8) is frictionless; i.e., p = 0. 

The seal deflection from the analysis results is computed by subtracting the seal gap (see Figure 11) at the end of 
the preload step from the seal gap at the end of the pressure step, or 


^ ^ pressure _ ^preload 


(i) 


The seal deflection is computed at two locations: directly in line with the center of the bolt hole and exactly 
between two bolt holes (i.e., a node on the plane of symmetry). Then the values from the two locations are 
averaged. 

The results for A z vs. pressure for the Baseline-ET model are compared to the average seal deflections from the 
ET test in Figure 12. As expected, the slope of the curve for the Baseline-ET model results is nearly linear. Thisis 
in contrast to the test results, which show nearly zero deflection up to 25 percent of the target pressure and then a 
linear response up to the target pressure. The reason for this bilinear response is hypothesized to be due to sticti on 
and load redistribution sometimes observed on the initial application of load to structural assemblies. Thetest report 
(Gillespie, 1988) indicates that the pressure was simply applied in 25-percent increments up to the target pressure 
and that deflection readings were taken (as discussed previously) at each increment. However, had the test 
procedure required applying only the first and second increments of pressure, then unloading the test specimen and 
re-zeroing the instrumentation, followed by applying the pressure increments to thefinal target pressure, it is 
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Normalized Seal Deflection 



Figure 11. Seal gap 5^ ressure due to the pressure load. 
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Figure 12. Comparison between test and baseline analysis results for seal deflection vs. pressure. 
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possible that a more linear deflection response would have been observed. However, this hypothesis cannot truly be 
verified unless a new test is performed. So while the slope of the predicted seal deflections does not match the test 
results, seal deflection values at the target pressure can be used to correlate the 3D analysis models to the ET test 
results. 

For the target pressure, the 3D ABAQUS analysis predicts a relative vertical deflection that is eleven percent 
different from the ET test result. Next, parametric studies are performed to determine the effects of three aspects of 
the FE modeling on the analysis results and establish better correlation to the test result. These parametric studies 
include 


• Effect of washer-beari ng-surf ace si ze 

• Effect of bolt-engagement BC 

• Effect of coefficient of friction 

Washer-beari ng-surface size. As shown in 

Figure 9(b), the spider constraint on the MHC 
connects the node that represents the bolt head to 
specified nodes on the outer surface of the MHC 
that represent the si ze of the washer beari ng surface. 
In the FE model, there are two rings of elements 
around the bolt hole. The outer ring has a diameter 
that is equal to the dimension of the washer’s outer 
diameter. In the Baseline-ET model, both rings of 
elements are used, as shown in Figure 13, 
representing the entire washer bearing down on the 
MHC when the bolt preload is applied. Two 
additional cases (see Figure 13) are also considered: 
(W1) the inner ring of elements is used for the 
spider constraint, and (W2) no rings of elements are 
used; i.e., the ring of nodes on the edge of the bolt 
holeisused for the spider constraint. 

Bolt-engagement BC. As shown in Figure 9(c), 
the spider constraint on the fixture connects the 
node that represents the bolt end to specified nodes 
on the inside surface of the threaded insert hole in 
the fixture thus representing the bolt-engagement 
BC. In the Baseline-ET model, the entire bottom 
surface of the hole is used, as shown in Figure 14. 
Two additional cases (see Figure 14) are also 
considered: (B1) the spider constraint at the bottom 
of the hole is maintained, and the ring of nodes at 
the mid-depth of the hole is used to define an 
additional spider constraint, and (B2) the spider 
constraints at the bottom and mid-depth of the hole 
are maintained, and the ring of nodes at the quarter 
depth of the hole is used to define an additional 
spider constraint. 



Figure 13. Cases considered to study the effects of 
washer-bearing-surface size. 


Figure 14. Cases considered to study the effects of bolt- 
engagement BC. 


Coefficient of friction. In the Baseline-ET model, a coefficient of friction of p = 0.1 is used. Two additional 
cases are also considered: (FI) p = 0.04, a general value for the coefficient of friction of Teflon against various 
materials, and (F2) p = 0.2. 

The results for seal deflection for the parametric studies are compared to the seal deflections from the ET test in 
Table2. lnTable2, A‘ est is the average seal deflection at the target pressure from the ET test, and all the deflections 
from the analysis predictions are normalized by this value. As expected, the seal deflections increase as the size of 
the washer bearing surface decreases. For the bolt-engagement BC, the additional spider constraints in cases B1 and 
B2 effectively shorten the bolt, and it is expected that the seal deflections would decrease. However, the results in 
Table 2 indicate that the seal deflection increases as the bolt is shortened. This result is not intuitive and warrants 
further study. Of the three parameters studied, the washer-bearing-surface size has the greatest effect on A z (which 
is consistent with results reported by Knight, et al. (2008)). The predicted value of A z for caseWI istwopercent 
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Table2. Seal deflections at the target pressure. 



ETTest Result 

ET BOSOR5 

ET ABAQUS Prediction 


Prediction 

Parametric Study 

Prediction 




Baseline- ET 

0.88 




W1 

0.98 




W2 

1.15 

A z 

1.0 

0.81 

B1 

0.91 

A‘ est 

B2 

0.93 




FI 

0.88 




F2 

0.89 




W1 +B1 

1.02 


different from the test result. The assumed bolt-engagement BC has a lesser effect on A ZI and the assumed 
coefficient of friction has a negligible effect on A z . 

One final case is considered that combi nes cases W1 and B1; i.e., one ring of elements on the MHC is used for 
the si ze of the washer beari ng surface, and two spi der constrai nts, one at the bottom and one at the mi d-depth of the 
threaded insert hole are used for the bolt-engagement BC. The seal deflection for this case isalso presented in Table 
2 and is less than two percent different from the test result. 

For comparison, the BOSOR5 analysis of seal deflection in the ET test was only within 19 percent of the test 
result (see Table 2). The ABAQUS analytical conf igurati on that combi nes cases W1 and B1 correlates most closely 
with, and conservatively over-predicts, the test result and hence will be used as the modeling practice for all 
subsequent analyses of the Ares-I US LH2 tank MHC seal joint. 

2. Ares-I US LH2 Tank MHC Seal Joint Analyses 

Through testing, the ET Project derived their own seal deflection allowable A e- for evaluating analysis results 
To minimize developmental testing, the US Project has decided to revert back to the original seal deflection 
specification. Therefore, in this paper, the seal deflections predicted for US are evaluated against the stricter design 
specification for seal deflection A,** The normalized seal deflections, (A z /A spec ), are reported. If (A Z /A^ )BC ) is less 
than unity, then the integrity of the seal is successfully maintained. If (A z /A spec ) is less than 0.5, then the seal 
deflection also satisfies the design specification with the one-sided interpretation. 

a. Finite Element Model 

The forward portion of the Ares-I LH2 tank is shown in Figure 15. The Naflex seal sits between the MHC and 
the dome-cover flange. The dome-cover flange is welded to the tank dome. The dome is connected to the tank 
barrel attheY-ring. The Y-ring is assumed to be rigid. 
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The MHC is shown in Figure 16. Most of the 
acreage of the MHC is a thin shell. To simplify the 
analyses, the thick regions at penetration bosses 
shown in the figure are ignored so that symmetry can 
be exploited. A wedge section of the MHC that 
encompasses a single bolt hole, or 3.913°, as shown in 
Figure 16, is used to construct the 3D FE model of the 
MHC. To construct the model of the Ares-I US LH2 
tank MHC seal joint, similar 3.913° portions of the 
seal, dome-cover flange, and dome are also used, and 
the modeling characteristics used for the ET test set- 
up anal yses are empl oyed. 

The 3D symmetric wedge FE model for the Ares-I 
US LH2tank MHC seal joint is presented in Figures 
17-21. The model for the enti re seal j oi nt assembl y i s 
shown in Figure 17. Note the locations of the Y-ring, 
the seal joint, and the tank axis. The model for the 
MHC is shown in Figure 18. The model for the seal 
is the same model used for the ET analyses and is 
shown in Figure 7. The model for the dome-cover 
flange is shown in Figure 19, and the model for the 
dome is shown in Figure 20. The dome-cover flange 
and dome are connected usi ng a mesh tie constrai nt to 
simulate the weld. A mesh tie constraint for shell elements couples all six degrees of freedom between nodes. 
A I ong the ti ed edge, the mesh on the dome has the same ref i nement as the mesh on the dome-cover f I ange, resul ti ng 
in a one-to-one nodal correspondence for the tie constraint. The entire model contains approximately 22000 
elements, including 21346 brick elements, 842 shell elements, and 10 beam elements 

Linear, elastic material properties are assumed for all of the materials The materials used for the current 
anal yses are provi ded i n Tabl e 3. 



Tank Axis 



Figure 17. 3D symmetric wedge finite element model of the Ares-I US LH2 tank MHC seal joint. 
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Solid Elements 



Shell Elements 


Shell-to-Solid Coupling 


Solid Elements 


Bolt Hole 


Shell-to-Solid Coupling 

Figure 18. 3D symmetric wedge finite element model of the Ares-I US LH2 tank MHC. 


Insert Hole 




Solid Elements 


Shell Elements 



Tie Constraint 
to Dome 

Figure 19. 3D symmetric wedge finite element model of the Ares-I US LH2 tank dome-cover flange. 



Figure 20. 3D symmetric wedge finite element model of the Ares-I US LH2 tank dome. 
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Table 3. Linear, elastic material 

properties for the Ares-I USLH2tank MHC seal joint anal' 

/ses. 

Part 

Material 

E, psi @ RT 

v @ RT 

MHC, Dome-Cover Flange, 
Dome 

Aluminum 2195 

10.7x10® 

0.33 

Seal 

Inconel 718 

30x10® 

0.3 

Bolt 

A286 Stainless Steel 

30x10® 

0.29 


The boundary conditions and applied loads are shown in Figure 21. The symmetric boundary conditions and 
contact definitions from the ET test set-up analyses are also used here. In addition, to simulate the rigid the Y -ring, 
a fixed boundary condition (i .e., u r = w e = u z = 9 r = 0 e = 0 2 = 0) is applied to the bottom of the dome. The bolts are 
assumed to be 180 ksi A286 fasteners, and the preload is taken to be 65% of the allowable yield stress of the bolt. 
The pressure load is applied along the entire inside surface of the assembly and around the pressure side of the seal 
legs as indicated in the Figure 21 insets. The pressure load is 1 .4 times the expected flight pressure to account for 
the joint-separation factor of safety requirement for US. Because there is a hole in the center of the MHC (see 
Figure 16) where a pressure system component is mounted, there is a thrust load associated with the pressure load. 
In this paper, the thrust load is assumed to be small because the diameter of the hole in the center of the MHC is 
small in comparison to the diameter of the MHC (D h0 ie/D MH c < 0.04). Therefore, in this paper, the thrust load is not 
included. 

The bolted connection is modeled with three spider constraints as established by the ET test set-up analyses. As 
such, for the bolt-to-MHC connection, one ring of elements around the bolt hole is used to represent the assumed 
bearing size of the washer even though the physical diameter of the washer is larger. In addition, for the bolt-to- 
dome-cover-flange connection, two spider constraints, one using the entire bottom surface of the insert hole and one 
usi ng the ri ng of nodes at the i nsert mi d-depth, are used to represent the bolt-engagement BC. 

The model described in this section will be referred to as the Baseline-US model. 



Fixed BC: 
u r =u 9 = u z = 0 

e r =e 0 = e z =o 


Tank Axis 


t 


Pressure Load 


Symmetry BC along cut planes: 
Solid Elements: u 0 = 0 
Shell Elements: 0 r = 0* = 0 


Figure 21. Boundary conditions and load applied to the 3D symmetric wedge finite element model of the 

Ares-I US LH2 tank MHC seal joint assembly. 
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b. Results 

The u z displacements for the Baseline-US model subjected to 1 .4 times the flight pressure are compared to the 
original unloaded configuration in Figure 22. (The FE discretization is excluded for clarity.) Because the dome is 
i ncl uded i n the model , the enti re seal j oi nt di spl aces, especi ally in the^ di recti on. A close-up of the di spl acements at 
the seal primary sealing surface is shown in Figure 23(b). The predicted average seal deflection A 2 (see Equation 1) 
for this model is (A z /A,^) = 0.39. This seal deflection is within the design specification with the one-sided 
interpretation. 


Seal Joint 



Figure 22. u z displacements for the Baseline-US model subjected to 1.4 times the flight pressure. 



Figure 23. Close-up of u z displacements at the seal primary sealing surface for the Baseline-US model 

subjected to 1.4 times the flight pressure. 
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The tank dome deforms axially and radially due to the pressure load. The M HC is free to move with the dome 
since, in the Baseline-US model, the MHC is restrained only by the symmetry boundary condition, the spider 
constraint that simulates the bolt, and the contact with the seal. In the assembled US, a pressurization system 
component would be bolted to the boss at the center of the MHC. Thi s component is not represented i n Basel i ne-US 
model. To represent the additional stiffness provided by this component, an additional radial BC is considered, as 
shown in Figure 24. The nodes along the inner diameter of the MHC are restrained in the r direction (i.e., u r = 0). 
The results for this case are identical to the results for the Baseline-US model. 

Radial restraint (u r = 0) 



Again, in the Baseline-US model, the MHC is restrained only through the spider constraint and the contact with 
the seal, and the dome deforms due to the pressure load. Because of the freedom of movement resulting from this 
configuration, the seal slides considerably from its original position, asshown in Figure23. This is in contrast to the 
ET test set-up model, in which the seal does not slide. For this reason, even though the results from the ET test set- 
up model indicate that the effects of coefficient of friction are negligible, the study is revisited here. Two cases are 
considered: p = 0.04, a general value for the coefficient of friction of Teflon against various materials, and p = 0.2. 
The results for the predicted average seal deflection for each of the cases considered are (A z /A spec ) = 0.39 and (A z 
/Aspec) = 0.38 for p = 0.04 and p = 0.2, respectively. These results are nearly identical to the results for the Baseline- 
US model. Theassumed coefficient of friction hasanegligibleeffect on predicted seal deflection. 

c. Cryogenic Temperatures 

Cryogenic temperatures are considered since differences in the coefficients of thermal expansion (CTEs) for the 
joint materialscan relieve the initial preload, resulting in greater seal deflection. 

The analysis is performed using the Baseline-US configuration and in three solution steps. In the first step, the 
bolt preload is applied. In the second step, the cryogenic temperature of -423°F is applied to the entire model. In 
the thi rd step, the pressure I oad i s appl i ed. 

An “at-temperature”-type analyses is performed. In at-temperature analyses, the material properties at the final 
temperature are used for the entire solution sequence (i.e., the material properties are not temperature-dependent). 
The materi al properti es at cryogeni c temperature used for thi s anal ysi s are provi ded i n Tabl e 4. 


Table 4. Linear, elastic material properties at cryogenic temperature. 


Part 

Material 

^cryo? pSi 

^crvo 

(Arvo, in/in/°F 

MHC, Dome-Cover 
Flange, Dome 

Aluminum 2195 

11.9x10® 

0.33 

8.8x10® 

Seal 

Inconel 718 

31.8x10® 

0.25 

5x10® 

Bolt 

A286 Stainless Steel 

30x10® 

0.29 

5.8x10® 


For the Baseline-US model, the bolt preload is based on the allowable yield stress of the bolt. For the cryogenic 
analysis, the load in the bolt decreases during the analysis due to the application of the temperature. The predicted 
average seal deflection is (A z /Aspec) = 0.28, which is less than the (A z /Aspec) due to the same pressure for the 
Baseline-US anal ysi sat room temperature (i.e., 0.39). 

Even though the preload decreases, the seal deflection decreases (in comparison to the Baseline-US model), 
which isan unexpected result. Subsequent analyses showed that this result is due to two factors. First, the increased 
stiffness of the materials at cryogenic temperature makes the individual seal joint components stiffer and more 
capable of resisting deformation. Second, the difference in the CTEs between the aluminum and Inconel causes the 
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MHC and dome-cover flange to compress against the seal primary sealing surface even further when the 
temperature is applied. Apparently, the effects of the reduction in bolt preload due to the cryogenic temperature are 
overcome by the response of the other components of the seal joint (the MHC, seal, and dome-cover flange) to the 
temperature. However, the result is still not intuitive and will be studied further. With the current results, the 
predicted seal deflection for the MHC seal joint assembly with cryogenic temperature is within the design 
specification with the one-sided interpretation. 


IV. Concluding Remarks 

Analyses of the Ares-I Upper Stage (US) liquid hydrogen (LH2) tank manhole cover seal joint were performed 
to assess the feasibility that the Naflex seal design used with the Space Shuttle External Tank (ET) manhole covers 
can successfully be used for the US manhole cover design. The analysis resultsfor the average seal deflections were 
compared to the US Project design specification limit for seal deflection. The US analyses were performed at 1 .4 
times the expected flight pressure to comply with the Constellation Program factor of safety requirement on joint 
separation. The seal deflections for all models considered are within the design specification; therefore, the US 
manholecover design using the Naflex seal should perform adequately. 

All of the analyses in this paper are three-dimensional symmetric wedge finite element analyses performed using 
the ABAQUS commercial finite element software. The analysis technique was validated by first modeling a 
heritage ET seal development leak test and correlating the analysis results for seal deflection to the test results. The 
modeling practices established using this ET model were used to construct the finite element models for the US 
analyses. 

Results from parametric studies performed with the ET models and the US models demonstrate the effects of 
several modeling assumptions on the average seal deflection. The assumed size of the washer bearing surface 
significantly affects the seal deflection, and the bolt-engagement boundary condition affects the seal deflection to a 
lesser extent. The assumed coefficient of friction for the contact definition has negligible effect on the seal 
deflection. 

The analyses described in this paper have established modeling practices that can be used for future analyses of 
the Naflex seals in other Ares-I US joint configurations. In these future analyses, the effects of preload and 
additional loading conditions on seal performance can be further evaluated. The use of commercial analysis 
software offers a great deal of flexibility both in modeling and interrogation of the results in comparison to the 
heritage axi symmetric modeling technique employed by the ET Project. 
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